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TECHNICAL PUBLICATION 


PROSPECTS FOR NUCLEAR ELECTRIC PROPULSION USING CLOSED-CYCLE 
MAGNETOHYDRODYNAMIC ENERGY CONVERSION 

1. INTRODUCTION 

From the most fundamental perspective, meaningful scientific and human exploration of deep space 
will be possible only if there is an abundant source of energy available for propulsion and onboard electri- 
cal power production. Furthermore, this energy must be stored in an extremely compact form if it is to be 
transported along with the spacecraft since any additional mass increment incurs a corresponding increase 
in energy requirement for orbital transfer. These fundamental considerations place a severe constraint on 
the minimum acceptable energy density needed for any deep space mission. 

To our knowledge, the only plausible means of providing the necessary energy densities, within the 
context of currently accepted physics, are beamed energy, nuclear energy (fission and fusion), and matter- 
antimatter annihilation. It should be noted that solar cells may be considered as the onboard conversion 
component of a beamed energy system with the caveat that the current specific mass for this technology is 
=30 kg/kW e and the perfected value is not expected to fall below 5 kg/kW e . The 1/r 2 decrease in solar 
radiation intensity also precludes robust missions beyond Mars. 

In principle, any of these energy sources would suffice from the standpoint of energy density alone. 
However, there are other attributes such as technical maturity, system complexity, energy conversion effi- 
ciency, and system specific mass, which stand to differentiate the various technology readiness levels. In 
terms of these broader system level considerations, nuclear fission space powerplants are the clear standout 
candidate for near-term realization. This train of thought is by no means unique and can be pursued to 
certain logical conclusions, as expounded most recently by Sackheim et al. 1 

Given this basic premise, the natural inclination has been toward nuclear-driven, closed-cycle gas 
turbine installations, since they represent the nearest evolutionary point beyond the existing technology 
base. Turbogenerator cycles, however, are only effective at low-heat rejection temperatures and as a result 
require large and massive space radiators. This performance characteristic derives from the fact that prac- 
tical material limits place severe restrictions on the maximum turbine blade operating temperature, even 
when the blades are actively cooled, and this, in turn, restricts the maximum cycle temperature. Despite 
this drawback, these systems have good prospects for approaching a powerplant specific mass below the 
10-kg/kW e mark. 

The major culprit affecting heat engine weight is the space radiator. Because the specific mass of 
the radiator is proportional to 1/T 4 , the radiator weight can greatly exceed the reactor and energy conver- 
sion component weights when forced to operate at a low-heat rejection temperature. Therefore, there is a 
strong incentive to break through the l-kg/kW e barrier by finding an alternative means of elevating the 



cycle operating temperature and thereby increasing the heat rejection temperature. The natural alternative 
for attaining this goal is magnetohydrodynamic (MHD) energy conversion. 

The primary reason for considering MHD is the ability to extract energy from a flowing gas at 
arbitrarily high temperatures. As such, an MHD converter can operate effectively at temperatures beyond 
the realm of actively cooled turbine blade technology. Using demonstrated solid core and fixed particle bed 
reactor technology, for example, it is possible to push peak operating temperatures to 2,500 K with 
minimum development risk, and this value could readily be raised to 3,000 K with further development of 
advanced high-temperature fissile fuels, particularly ternary carbide elements. Even operating tempera- 
tures beyond 3,000 K are conceivable with the introduction of gas core reactors in which the fissile fuel is 
allowed to circulate throughout the power conversion system. The useful upper limit for gas core reactors 
with MHD energy conversion is estimated to be as high as 8,000-10,000 K, at which point the rapid 
increase in free electron radiation makes regenerative cooling extremely difficult as well as causing the gas 
cooling timescale to fall — most likely below the channel residence timescale. 

The most common configurations for MHD generators are the (a) linear channel and the (b) radial 
disk, as illustrated in figure 1. In either case, an external magnetic field, B . is imposed perpendicular to the 
flow path and the motion of the electrically conductive gas interacts with the magnetic field to develop a 
Faraday electromotive force (emf) orthogonal to both the flow velocity and the magnetic field and a Hall 
emf directed along the flow path. The value of the Hall emf relative to the Faraday emf depends upon the 
value of the Hall parameter, (i When /I is on the order of unity or greater, the Hall configured disk geometry 
is used predominately since it avoids the complicated electrode segmentation requirements associated with 
the design of linear channels. 


(a) 


Segmented Anode 




Figure 1. Schematics of (a) linear and (b) disk MHD generator configurations. 
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The induced Faraday and Hall emf ’s result in a current flow in the gas, which is transferred to an 
external load through contact electrodes in the channel. The interaction of the current with the magnetic 
field produces a retarding force on the gas, thereby forcing it to do push work as it moves through the 
generator channel. Therefore, there is a pressure drop across the MHD converter in direct analogy to that 
experienced by a turbine stage. The difference being that the stationary and lightly stressed MHD channel 
walls can be cooled to temperatures much lower than the core gas temperature in the same manner that the 
chamber walls of a chemical rocket engine are regeneratively cooled. In this case, it is the heat transfer rate, 
not the core gas temperature, that determines engineering feasibility. 

Nuclear MHD space power systems, as summarized above, were first conceived and conceptually 
developed by R J. Rosa in the 1960’s. 2-8 Since then, the idea has been periodically revived and amplified 
in relation to various perceived space power needs. 9-15 Yet, this technology has never progressed beyond 
studies and low-priority research efforts, despite its many positive attributes. And it remains to this day, a 
potentially high payoff concept for which the perceived risk far exceeds the actual development risk. 

This Technical Publication revives the concept once more and amplifies on the scenario for a nuclear 
MHD space powerplant coupled with high-power electric thrusters. The emphasis here is to demonstrate 
that these systems could immensely broaden deep space mission capability while adhering to a technology 
development path that is entirely within established bounds of engineering feasibility. These conclusions 
are obtained through a preliminary but thorough system analysis in which the various component charac- 
teristics are estimated using established data or, when an empirical database is unavailable, using widely 
accepted extrapolations from current state-of-the-art technology. Additional technological development 
opportunities are also examined in terms of potential payoff versus development risk. 
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2. SYSTEM ANALYSIS 


Generally speaking, high-power density fission reactors are adaptable to both Brayton and Rankine 
energy conversion cycles. For example, Rankine cycles, by taking advantage of the gas-to-liquid phase 
transition, provide superior performance for a fixed overall temperature differential, and they can com- 
pletely avoid the limitations of compressor turbomachinery. Unfortunately, the condensing vapors tend to 
be highly corrosive, raising significant design and reliability concerns. Brayton cycles, on the other hand, 
are simpler, exhibit only slightly inferior performance, and promise increased reliability in that they are 
more suited for utilization of an inert gas working fluid. 

Although previous studies have established the fact that Rankine MHD cycles can provide substan- 
tially lower system specific mass, 2-8 ’ 11 ’ 12 the authors take the point of view that Brayton MHD cycles are 
more directly compatible with solid-core, gas-cooled reactors and can more effectively leverage the exist- 
ing technology base. Therefore, they appear to offer a lower risk development path in the near term while 
maintaining a significant payoff in comparison to turbogenerator systems. As such, the Brayton MHD 
cycle with regenerative heating and intercooled multistage compression forms the focus of attention for 
further investigation. 

The essential components of a nuclear electric propelled spacecraft based on a closed-loop MHD 
Brayton cycle are shown in figure 2. From the perspective of mass distribution, it is useful to lump the 
spacecraft subsystems into three categories: the payload, the propellant, and the thrust producing system. 
The major thrust of this work was to analyze the specific mass characteristics of the powerplant and the 
integrated propulsion system and to make a status/risk assessment of the various enabling technologies. 
Since the powerplant is the dominant contributor to overall spacecraft weight, there is intense interest in 
obtaining accurate estimates of powerplant specific mass characteristics as a means of delineating 
technological capabilities and limitations. Once the specific mass characteristics are quantitatively defined 
to a satisfactory degree of accuracy, it becomes possible to comparatively assess various technological 
approaches and to identify the major technology development challenges and associated risks. 
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Figure 2. Essential components of a nuclear electrical propelled spacecraft 
based on a closed-loop MHD Brayton cycle powerplant. 
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For reference purposes, a schematic of a closed-loop nuclear MHD Brayton cycle powerplant coupled 
to high-power electric propulsion is provided in figure 3. The temperature-entropy diagram for the closed- 
loop Brayton power cycle is also shown in figure 4. The thermodynamic cycle is now analyzed assuming 
fixed maximum (T max ) and minimum (T min ) operating temperatures. 


MHD Generator 


Auxiliary 

Power 





0 


rad 


Figure 3. Schematic of a closed-loop nuclear MHD Brayton cycle powerplant 
with regenerative heating and intercooled multistage compression. 
Installation is shown coupled to high-power electric propulsion. 
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Figure 4. Temperature-entropy diagram for a closed-loop Brayton cycle 
powerplant with regenerative heating and intercooled multistage 
compression. 


2.1 Thermodynamic Cycle Analysis 


In the special case where all net power production is used for propulsion, the thermal efficiency 
(?] th ) of the powerplant is equivalent to the ratio of thruster electric power (W t hru S te r ) to reactor thermal 
power (Q reactor ): 

„ ^thruster VV g en ~ Wcom P ( r max ~ ' T 9 ) ~ N c ( 7 2 ~ 7 min ) 

• th S-\ S-\ ry rj i ? ^ / 

'■^reactor ^reactor 1 max 1 7 

where IF gen is the electrical power output of the generator, W comp is the power consumption of the com- 
pressor bank, and N c is the number of compressor stages. In addition, it is assumed T min = 7^ = T 3 = T 5 . . . 
is the effective radiator temperature. 

In analyzing the MHD generator, it is useful to note that any MHD device can be described in terms 
of two fundamental process parameters: an enthalpy extraction/addition ratio (r\ N ) and an isentropic effi- 
ciency ( i] s ) which quantifies the degree of departure from an isentropic process. 

The enthalpy extraction parameter for the generator is defined as the ratio of the change in enthalpy 
A /t to the entrance enthalpy h ent : 

_ A^gen _h%- h 9 _ T max - T 9 (2) 

TIn ~ ~h _ — / = — T 

"cnl "8 max 
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or 


^9 (l max • (3) 

For an MHD generator (enthalpy extraction), the isentropic efficiency rj s is defined as the ratio of 
the actual change in enthalpy to that for an ideal isentropic process Ah gen s : 

_ Ah gen hg ~ hg _ T max - Tg (4) 

Ahg Cns hg — hg s T max — Tg s 

or 

r 9 =| l ~ r ls,g 

Using the isentropic process relationship for an ideal gas, it is possible to write equation (5) in the form: 


'9, 


L max 


L max 


(5) 




y- 1 " 



II 

1 ^S,g 

\-n g '/ 

' 7inax ’ 

(6) 


where n g = Pg / Pg is the pressure ratio across the generator and yis the specific heat ratio. Combining 
equations (3) and (6) yields a fundamental relationship between r\ N and ii s : 


Vn=V . 


s,g 


y-i 

l ~ n g 7 


(7) 


The N c compressor stages are assumed to operate with the same temperature rise AT C = T c - T min 
such that T c - T 2 = T 4 = T () ... Furtheimore, in order to compensate for the net pressure drop in the system, 
the pressure ratio across the MHD generator is defined as jv g = Gn^ c such that each compressor stage has 
a common compression ratio (n c )\ 


n 


C 




1 /N c 


= P2 
Pi 


PA_ = P6 
P3 P5 


( 8 ) 


The parameter G represents a pressure loss ratio due to frictional losses (0<G<1). Thus, for each 
compressor stage (enthalpy addition), the isentropic efficiency ( /) s c ) is defined as 
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'^corrip h<: ^rniri -Tc ^min 


(9) 


which can be reformulated in terms of T c : 
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( 10 ) 


By using the isentropic process relationship for an ideal gas and using equation (8) to eliminate^., 
it is possible to write equation (10) in the form: 


t c =t 2 =t 4 = t 6 = 
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(ii) 


Regenerative heating is employed in the cycle to increase thermal efficiency and reduce the amount 
of waste heat that must be rejected by the radiator. The regenerator effectiveness is defined as 
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<7 max re g en ^9 ~ ^6 ^9 — ^6 
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or, upon substitution of equations (6), (7), and (11), 
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A power balance on the regenerator requires that 


h 9 ~ %0 _ t 9 ~T\() _ 1 (I 5 ) 

h l -1% Tj- T 6 
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or 

T m=T 9 -( r 7 -T 6 ) . (16) 

Thus, equations (13) and (16) may be combined to obtain the form : 


*10 — (l — £ regen ) *9 _ £ regen*6 • ( 17 ) 

Then, substitution for T 9 and T () using equations (6), (7), and (1 1) gives an expression for T ]() in teims of the 
cycle operating temperatures and process parameters: 
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The vaiious theimodynamic states throughout the cycle have now been expressed in teims of the 
reactor and radiator operating temperatures and the vaiious component process parameters. This informa- 
tion may further be used to deduce the power per unit mass flow per unit specific heat for each system 
component. The results are summarized as follows: 
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It is also possible to frame the cycle thermal efficiency in similar terms by substituting equations 
(19)— (21) into equation (1): 
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2.2 Specific Mass Analysis 

Thermal efficiency is not the central consideration for space installations. Rather, the entire system 
must be examined in terms of its overall mass relative to its performance. This requires a careful evaluation 
of all subcomponent mass characteristics in order to obtain an estimate of the aggregate system mass. 

The aggregate specific mass of the nuclear MHD powerplant is defined as 


aq . + awl + awl + aq + aq 

^ reactor 'gen icomt) 1 repen 1 
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plant 


comp 


i regen 


I rad 


(25) 


Wgen W’comp 


The specific mass of the combined powerplant and electric thruster, which is assumed here to 
include all necessary power conditioning equipment, is therefore given by 


^system ^plant ^thrust 


(26) 


It is also common practice to define the specific power as (j) = 1/a. Specification of the subcomponent 
specific mass characteristics now follows. 
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2.2.1 Fission Reactor 


Ample experience indicates that solid core and fixed particle bed reactors can provide peak tem- 
peratures in the 2,500-3,000 K range. The solid-core nuclear engine for rocket vehicle application (NERVA) 
engine, for instance, had a design operating temperature of 2,500 K, and a low-risk development path was 
identified for increasing this limit to 3,000 K through the development of high-temperature carbide fuels. 16 
As such, this technology base can be used to estimate anticipated reactor specific mass. Using the NERVA 
technology as an empirical basis, for instance, Holman has developed a 350 MW th reactor design having a 
mass of 1,785 kg. This yields a specific power of 0.2 MW th /kg. 

Within reasonable bounds, the reactor mass can be assumed independent of reactor power level, 
and the authors conservatively estimate a reactor mass of 3,000 kg, which includes margin for shielding. 
Thus, 

_ ^reactor _ 3,000 kg 

^reactor ~ r) ~ r\ • A ' ) 

^reactor k'rcaclor 


2.2.2 Magnetohydrodynamic Generator 

There are strong incentives to operate MHD devices at very high magnetic flux densities. For 
example, the power density of an ideal MHD generator varies as ou 2 B 2 , where o is the bulk electrical 
conductivity, u is the velocity, and B is the magnetic field strength. Since a and u have practical limitations, 
it is therefore desirable to have B as large as practically feasible. Furthermore, Hall -configured disk genera- 
tors can exhibit high isentropic efficiency only when the Hall parameter is large (j8 >10). 17 Again, high 
magnetic flux densities are required since /3 °c B/p, where p is the static pressure in the device. 

Because of this basic need for large values of B, the weight of the generator as a whole is dominated 
almost exclusively by the weight of the magnet. The magnet weight, in turn, is determined mainly by the 
weight of the confinement structure and the weight of the coils. Under these assumptions, the generator 
specific mass is defined by 


(Zg en ^magnet ^struc ®coil , ( 28 ) 

where a ina;2nct is the specific mass of the magnet, a stmc is the specific mass of the structure, and a coil is the 
specific mass of the coil windings. 

It is not practical to entertain the utilization of dissipative magnets in space installations, mainly 
due to the excessive weight penalties associated with waste heat rejection. Therefore, only supergenic 
(i.e., superconducting-cryogenic type coils) are considered in which a high field superconductor 
(e.g., M^/Sn or Nb/Ti) is embedded in a metal base material (i.e., copper or aluminum) as a means of 
providing quench stabilization. It is also assumed that flux pumping techniques will be utilized in order to 
minimize the auxiliary equipment needed for magnet power-up. 
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2.2.3 Magnet Confinement Structure 

The basic challenge of large-volume, high-field magnet design can be readily summarized. 
Because magnetic flux density falls off with distance from the coil as 1/r 2 , very large currents are necessary 
to fill the working volume, and this, in turn, leads to the exertion of extremely large Lorentz forces on the 
coil. Therefore, a large confinement structure is needed to support those forces that the conducting coils 
cannot withstand themselves. 

The confinement requirements may be fundamentally expressed in terms of the stored magnetic 
field energy, W m : 



dV = 


^—v ’ 
2 Po 


(29) 


where p () is the magnetic permeability and V is the enclosed volume. The Virial theorem may then be used 
as an estimate of the minimal mass due to structural requirements (i.e., ideal hoop tension to contain the 
stored energy): 


Pi 

s t 


m stiuc — — 


St 2^0 


-V 


(30) 


where p is the material density and s t is the material working stress. Some typical values for W m /m = s t /p 
are summarized in table 1. 


Table 1. Virial theorem requirements. 


Material 

W m /m=s t /p( kJ/kg) 

Fiber-reinforced composites 

10-50 

Stainless steel (304LN) 

44 

Aluminum (221 9T851) 

107 

Titanium 

309 

Beryllium-Copper 

580 


At this point, it is convenient to introduce the generator power density, P gen = W„ cn /V, into 
equation (30) to arrive at a working specific mass estimate for the confinement structure: 


^struc 


w struc 

w 

,r gen 


P £ 2 /2pQ . 

s t Pgen 


(31) 
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2.2.4 Saddle Coil Winding — Linear Generator 


Consider the coil requirements for a linear generator channel. Based on conservative design expe- 
rience, assume a minimal channel aspect ratio of L/D = 10. This ratio is recommended by Rosa in order to 
minimize nonideal losses associated with large surface-to-volume effects. 18 In this case, the working 
volume of the generator can be expressed as 


V, 


W, 


gen 


gen 

o 

gen 


*d 2 l 

4 


5n ^ 
— D 5 
2 


which yields a nominal generator diameter of 


(32) 


D = 


2 W, 


gen 


5 n P. 


gen J 


(33) 


A uniform magnetic field can be produced within the generator volume if the saddle coil windings 
approximate a circular crescent cross section of diameter D, as illustrated in figure 5. The total cross- 
sectional area (A c ) of the winding on each side is 


A c ~ 2ca = 2c\ y + ^ \ = c(c + D) = cD\ 


D 


(34) 



Figure 5. Illustration of a circular crescent coil cross section for a linear channel 
of diameter D. Note that 2a = D + c is a geometric constraint. 


In the limit c/D —> 0, A c ~ cD and a uniform constant magnetic field (B c ) is produced which is 
defined by 


_ A) _ _ Ppjc^c _ Ppjc c , 

2 4r 2D 2 


(35) 


14 



where i is the total current, j c is the current density, and B () = q () //2r is the field produced in a loop of 
diameter D = 2 r. The value for j c can be as high as 10 9 A/m 2 for Nb 3 /Sn-based superconductors. 

Using equations (33) and (35) to eliminate D and c in equation (34) gives an expression for the 
cross-sectional area in terms of the material properties, the generator power characteristics, and the 
required magnetic flux density: 


, 2 B c 

A: = “ 

BOJc 


2 B<. 
Bo Jc 


■ + 


2 W, 


gen 


y57t -fgen J 


The mass of the magnet coils may now be deduced from the relation: 

tn c — p c V c — 2 PcA^ ~ 20 p c A c D 
and the specific mass of the coil windings takes the form: 


°xoil 


m c _ 20p c A c D 


W W 

rr gen vv gen 


(36) 


(37) 


(38) 


where D is defined by equation (33). 

2.2.5 Helmholtz Coil Winding — Disk Generator 


Consider the coil requirements for an MHD disk generator of radius r and mean channel height h. 
Based on conservative design experience, assume a minimal channel aspect ratio of r/h =10 in analogy to 
the length-to-diameter ratio assumed for the linear configuration. Again, this value is recommended by 
Rosa in order to minimize nonideal losses associated with large surface-to-volume effects. 18 Thus, D/h 
= 20, and the working volume of the generator can be expressed as 


V. 


W, 


gen 


gen 


gen 


7t , 7T Q 

— D z h = — D' 
4 80 


(39) 


This yields a nominal disk diameter of 


D = 


SOW. 


gen 


n P. 


gen J 


(40) 
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Helmholtz coil windings of diameter D, as illustrated in figure 6, can then be used to produce a 
uniform magnetic field B H within the disk volume: 

R B 0 M /'p/r-'V . (41) 

" V2 2-n r -J2D 

This expression can be combined with equation (40) to arrive at an expression for the cross- 
sectional area of the winding in terms of the material properties, the generator power density, and the 
required magnetic flux density: 




Ac 

2 



Figure 6. Illustration of a Helmholtz coil winding for a disk generator of diameter D. 
Note that the coils are separated by a distance r = D/2. 




V2 B 


H 


BoJc 


80 W. 


gen 


V 71 ^gen ) 


The mass of the magnet coils may now be deduced from the relation: 


m c PcV c Pc nD A: 


427ip c Bj]D~ 

PoJc 


and the specific mass of the coil windings takes the form: 


°xoil 


m c _ 42np c B H D~ 


W, 


gen 


PoJcW, 


gen 


(42) 


(43) 


(44) 


where D is defined by equation (40). 
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2.2.6 Turbocompressors 


The weight characteristics of the turbocompressors and drive can be accurately estimated from a 
broad body of technical experience. Following Rosa, it is assumed that the specific mass of the compressor 
group is given by 


^corrip 


^» = 2xlO- 5 kg/W . 

^comp 


(45) 


2.2.7 Regenerator 

The regenerator is assumed to be a typical compact shell and tube heat exchanger for which the 
specific mass attributes can be estimated accurately. A conservative value for the mass per unit area of this 
type of heat exchanger is 


Aege„=y S& = lkg/m 2 . (46) 

A'egen 

Assuming a constant value for the overall heat transfer coefficient, t/ regeri , it is possible to express 
the total heat transfer rate in the form: 

2regen ~ ^ regen A^ lMD A regen > (47) 

where A!T LMD is the log-mean temperature difference for the heat exchanger and is defined as 


and 


In the limit I AT b /A T a | — > 1 , 


A7i 


A T b - A T a 


LMD - / AT \ 

A lb 


In 


v A T a J 


AT a 7jiot,in ^cold.out 


A Tb ^hot.out -^cold.in • 


limi 


A T b 
A T n 


AT \ 


LMD 


A Tb + A T a 


-»1 


(48) 


(49) 

(50) 


(51) 
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Note that the temperatures used to calculate AT l MD are obtained from the thermodynamic cycle analysis 
based on an assumed value for the regenerator effectiveness. 


If equations (46) and (47) are now combined to eliminate A c , the specific mass of the regenerator 
takes the form: 


a 


regen 


m regen _ Aegen 
Qrcgcn ^regen^LMD 


(52) 


2.2.8 Radiator 

Space radiators tend to dominate the weight of multimegawatt heat engine installations, and nuclear- 
electric propulsion systems can quickly become impractical due to the impact on overall system specific 
mass. Issues of reliability and suivivability also escalate with increasing size. The important design param- 
eter for space radiators is the mass per unit area defined as 

Aad=T^ = l — 0.2kg/m 2 . (53) 

Aad 

The upper range for /3 rad represents projected estimates for conventional radiator design whereas the lower 
range for /3 rad corresponds to projections for liquid droplet radiator concepts. 

The radiative power is given by 


Grad £C7AradT m in , (54) 

where eis the emissivity and a is the Stefan-Boltzmann constant. If equations (53) and (54) are combined 
to eliminate A rad , the specific mass of the radiator takes the form: 


(Z ra d 


,n rad 

Grad 


Aad 

£OTl. 


(55) 
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2.2.9 Electric Thruster 


Many nuclear-electric deep space transport missions exhibit optimal results for an idealized thruster 
having a specific impulse in the range of 2,500 to 1,500 sec and a thrust in the range of 10 to 1,000 N. 
For a thruster efficiency Tj t = 50%, the fundamental equation of electric propulsion would imply a nominal 
power range of 


W e =F^-^~ 0.25-15 MW e , (56) 

where F is the thrust, I is the specific impulse, and g 0 is the Earth’s gravitational constant. 

In reality, efficient, high-thrust electric propulsion has yet to be demonstrated, and the prevalent 
opinion among specialists is that electric thrusters can only be effective at low thrust levels. Nevertheless, 
there is some hope that certain electric thruster devices may be scaleable to multimegawatt power levels 
and yet remain compact and reliable. The variable specific impulse magneto plasma rocket (VASIMR) 
concept even holds the promise of constant power throttling which allows both thrust and specific impulse 
to be continuously varied. 19 Extensive calculations by Slavin et al. also indicate that an MHD high- 
temperature current- (T-) layer type of engine using neon propellant could potentially produce 1,000 N of 
thrust at a specific impulse of 2,000 sec with an electrical efficiency of more than 90%. 15 

Some representative electric thruster characteristics, including conservative estimates for specific 
mass, are summarized in table 2. Note that these are only approximate figures of merit based on a broad 
body of technical experience. Here, the specific mass is defined in terms of the plant electrical power 
delivered to the thruster: 


_ ^thrust _ w thrust 

t " n,SI - WVus, “W'gen-W'comp 


(57) 


Table 2. Electric thruster characteristics.* 


Device 

■8" 

'c/T 

CD 

O 

>h 

a (kg/kW e ) 

Ion (Kr) 

>5,000 

0.8 

1.0 

MPD (Li) 

4,000-8,000 

0.5 

0.5 

MPD (H 2 ) 

>8,000 

0.5 

0.5 

VASIMR (H 2 ) 

3,000-30,000 

0.5 

0.2-1 .0 


‘Approximated values 
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3. SYSTEM ASSESSMENT 


Based on the preceding development, representative system analysis calculations were performed 
as a means of assessing system attributes. First, thermodynamic cycle characteristics were investigated 
over a practical range of T min and 7C ; i.e., r] N , parameter space. The assumed system characteristics, based 
on established bounds of engineering capability, are as follows: Helium working fluid, T max = 2,500 K, 
£ regen = 0-9, V s g = 70%, and i] s c = 87%. All calculations were performed assuming no frictional pressure 
losses (i.e., G= 1). 


The resulting thenual efficiency contour map for the MHD Brayton cycle is shown in figure 7. Note 
that the peak cycle efficiency moves to larger pressure ratios as the minimum cycle temperature falls. The 
highest achievable thermal efficiency approaches 70% as the heat rejection temperature approaches 
200 K. For space applications, however, thenual efficiency is secondary in importance to system specific 
mass. 


MHD Brayton Cycle Thermal Efficiency 
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Figure 7. Thenual efficiency of an MHD generator Brayton cycle using helium as the working fluid. 


For demonstration purposes, the authors assumed a reactor power of W g&n = 100 MW th and 
explored the specific mass characteristics over a similar range of 7j r|in and (i.e., rj N ), parameter space. 
Two sets of subsystem technology assumptions were investigated. These assumptions were based on 
conservative estimates of near-term and advanced subsystem technologies, as summarized in table 3. The 
major anticipated technology advances are as follows: An increase in s t /p through utilization of beryllium- 
copper material for the magnet structure, a decrease in p c through utilization of aluminum as the base 
material for the supergenic coil, major reduction in radiator specific mass through the introduction of 
liquid droplet technology, and slight evolutionary improvements in thruster efficiency and specific mass. 
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Table 3. Subsystem technology assumptions. 


Parameter 

Near Term 

Advanced 

^reactor (k9) 

3,000 

3,000 

P gen (MW e /m 3 ) 

500 

500 

B (tesla) 

8 

8 

a/s f (kJ/kg) 

309 

580 

Pc (kg/m 3 ) 

1 xIO 4 

3 x 10 3 

k (A/m 2 ) 

1 xIO 9 

1 xIO 9 

U regen (W/m 2 K) 

500 

500 

Pregen (kg/m 2 ) 

1 

1 

e rad 

0.9 

0.9 

/l rad (kg/m 2 ) 

1 

0.2 

Othrust (°/°) 

50 

60 

a thmst (kg/kw e ) 

0.5 

0.4 


The resulting contour maps for the disk MHD powerplant specific mass are given in figure 8. For 
the near-term technology assumptions, it is found that the powerplant specific mass can be less than 
0.4 kg/kW e at a heat rejection temperature of =700 K and an enthalpy extraction ratio of =20%. For the 
advanced technology assumptions, the powerplant specific mass approaches 0.2 kg/kW e at a heat rejection 
temperature of =500 K and an enthalpy extraction ratio of =30%. Results for the linear MHD channel 
configuration display virtually identical quantitative behavior and are not shown. 
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Figure 8. Specific mass characteristics of a disk MHD Brayton cycle powerplant 
using helium as the working fluid. 
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In order to keep the powerplant’s circulating heat power comparable in magnitude to the generated 
electrical power, it is desirable to utilize the highest practical enthalpy extraction ratio consistent with 
specific mass optimization desires. Accepting a generator pressure ratio of n = 8 as an optimal value yields 
an enthalpy extraction ratio design value of t] N = 40% when 77 = 70%. The powerplant and system 

specific mass characteristics for this particular case are indicated in figure 9, assuming a target radiator 
temperature of 500 K. Note that the propulsion system specific mass is at the 1 kg/kW e mark or less using 
very conservative design parameters. Also observe that the thermal efficiency now assumes a respectable 
value of rj th = 40%. 



200 300 400 500 600 700 800 


7min(K) 



200 300 400 500 600 700 800 


7min (K) 

Figure 9. Disk MHD powerplant specific mass attributes for i] N = 40% 
and rj s = 70%. The working fluid is helium. 
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The predicted scaling of MHD powerplant and propulsion system specific mass with net electrical 
power generation, under the assumption that rj N = 40%, is shown in figure 10 for peak reactor 
temperatures of 2,500 and 3,000 K. These estimates turn out to be in substantial agreement with the scaling 
results obtained by Rosa and Myrabo. 12 The scaling relationship developed by Mathematical Sciences 
Northwest, Inc. (MSNW) for closed-loop Brayton cycle turbogenerator systems is also shown for 
comparison. 20 Although the MSNW correlation indicates an economy of scale with increasing power, it is 
clear that the MHD Brayton cycle promises much greater payoff for a similar level of development risk. 




Figure 10. Scaling of MHD powerplant and propulsion system specific 
mass with net electrical power. 
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4. MAGNETOHYDRODYNAMIC TECHNOLOGY ISSUES 


The success of closed-cycle nuclear MHD powerplants in space will depend on the ability to design 
reliable MHD generators with high- efficiency parameters. The previously discussed system analysis 
results, for example, depend on assumed values of 40% and 70% for the enthalpy extraction ratio and 
isentropic efficiency, respectively. In fact, powerplant specific mass is extremely sensitive to decreases in 
isentropic efficiency, and it is possible to assert that these types of installations will be competitive only if 
7] s g > 65%. 15 

A critical review of the technology reveals that the greatest hope for attaining these performance 
goals lies with nonequilibrium ionization Hall generators. Years of technology development effort in this 
arena, mainly by Japanese researchers, have yielded demonstrated efficiency parameters of i] N = 30% and 
T) s = 50%. Clearly, this level of performance is insufficient for space installations, but these results are 
based on small-scale experimental facilities, which suffer severe surface-to-volume penalties. Based on 
detailed numerical simulations and known scaling laws, these same researchers are confident that the 
efficiency parameters can be raised to the necessary level at practical generator scales. 

In this section, the conceptual design of a high-performance, seeded nonequilibrium disk generator 
is presented. Opportunities for exploiting neutron-induced ionization mechanisms are also discussed as a 
potential means of increasing electrical conductivity and enhancing performance and reliability. 

4.1 Conceptual Design of the Nonequilibrium Disk Generator 

Experience has shown that the successful design of a seeded nonequilbrium disk generator requires 
the simultaneous satisfaction of many conflicting restrictions. In addition to keeping the plasma stable and 
the isentropic efficiency high, there are practical limits on outlet Mach number, maximum wall divergence 
angle, maximum electric field strength, maximum electrode current density, and maximum effective elec- 
trode width. 

All of these considerations are important from an engineering design perspective. For example, the 
outlet Mach number should be only slightly greater than unity in order to obtain good diffuser efficiency, 
the maximum wall divergence angle should be maintained less than 10° in order to prevent boundary layer- 
separation, the maximum electric field strength and current density depend on available insulator and 
electrode materials, and the effective electrode width is always less than the actual electrode width due to 
current concentration. 

As previously noted, Hall-configured disk generators can exhibit high isentropic efficiency only 
when the Hall parameter is large 1/3 >10), implying the need for high magnetic flux density. 17 However, it 
is just as important to note that nonequilibrium MHD plasmas are subject to ionization instability. This 
instability can create inhomogeneous gas dynamic and electrical properties and lead to severe performance 
degradation. 21 Past research has shown that this can be avoided in a helium plasma with fully ionized seed 
if the temperature is maintained in the 4,000 K < T e < 7,000 K range. 22 
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To estimate the conditions necessary for satisfying this requirement, it is possible to utilize a sim- 
plified expression developed by Okuno et al. for a design parameter F defined as the ratio of the radial 
electric field strength to the magnetic flux density: 23 


F = 



R 


3t „-F2M1 To 


1 + 


7-1 


(58) 


where E r is the radial electric field, R is the gas constant, T e is the electron temperature, M is the Mach 
number, and T () is the stagnation temperature. 

A plot of this parameter is shown in figure 1 1 for typical values of M and T () . From this result, it is 
possible to deduce that the electric field strength per unit magnetic flux density must be in the range of 
4,000 V/m-T < F < 5,000 V/m-T for the plasma to remain stable. Thus, given a permissible range on 
stagnation temperature and electric field strength, specification of B and E r is sufficient to define plasma 
stability. 


Design Space for Cesium-Seeded Helium 



Te( K) 


Figure 11. Variation in design parameter F with M and T e 


Inui, Ishikawa, and Umoto have used this basic condition to develop a conceptual design methodol- 
ogy, which avoids the need for trial and error search procedures. 17 Using this methodology, they have 
developed a conceptual design for a 100-MW th disk generator consistent with space power installation 
requirements. Their assumed design restrictions, based on practical engineering considerations, are 
summarized in table 4. 
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Table 4. MHD disk design restrictions. 


Parameter 

Value 

^reactor (MWth) 

100 

To (K) 

2,200 

P 0 (atm) 

3 

/in (m) 

0.3 

%(%) 

40 

B (T) 

8 

A^exit 

1.2 

®max (°) 

10 

Er,max (kV/m) 

35 

Anax (kA/m 2 ) 

35 

Effective Electrode Width (m) 

0.025 


Their resulting design for disk lofting is shown in figure 12 where the outlet radius is 70 cm, the 
inlet height is 5.8 cm, and the exit height is 8.5 cm. For detailed variation in gas dynamic and electrical 
properties throughout the channel, the reader is directed to the original reference. 17 The major points to 
note are that the plasma is stabilized at an electron temperature of T e ~ 5,000 K and the i sen tropic efficiency 
is predicted to be rj s ~ 83%. 
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Figure 12. Design lofting for a 100-MW th MHD disk generator 
using cesium- seeded helium working fluid. 


The purpose of this conceptual design exercise is to demonstrate that MHD disk generators can be 
developed which exhibit the necessary performance criteria while simultaneously satisfying an array of 
practical engineering design restrictions. Unlike many advanced concepts, realization of nuclear MHD 
space power does not depend on any assumptions of idealized performance. Certainly, there are major 
challenges to be faced in developing this technology, but unlike many advanced concepts, one may proceed 
on a sound basis of realism in terms of both performance estimates and engineering feasibility. 
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4.2 Neutron-Induced Ionization Enhancement 


Historically, closed-cycle nonequilibrium MHD research has focused on cesium-seeded inert gas 
working fluids. However, nuclear reactor thermal sources present opportunities for nuclear-induced 
ionization processes, which may greatly enhance electrical conductivity and eliminate the need to circulate 
condensable seed material. 

In this ionization process, a working medium is either composed of or seeded with an isotope 
having a large neutron interaction cross section. The interaction of the neutron and the isotope results in the 
production of charged particles with kinetic energy, large compared to the ionization potential of the sur- 
rounding medium. These energetic charged particles then excite and ionize the surrounding molecules via 
collisions. Secondary ionization processes may also occur due to collisions between the liberated primary 
electrons and surrounding molecules. Some examples of interesting interactions with thermal neutrons are 
listed in table 5. 


Table 5. Large cross-section neutron interactions. 


Interaction 

Cross Section 
at 0.025 eV (barns) 

Kinetic Energy 
Release (MeV) 

3 He(n,p) 3 H 

5,300 

0.76 

10 B(n,a) 7 Li 

3,840 

2.30 

6 U(n,a) 3 H 

940 

4.78 

235 U(n ,ff)ff* 

580 

=170 


*ff denotes fission fragment 


The kinetic energy released in nuclear interactions such as these is at least four orders of magnitude 
larger than the ionization energy of any molecule. As a result, each nuclear interaction is capable of causing 
thousands of ionization events, assuming all of the kinetic energy is transferred to the surrounding 
medium. Two mechanisms that can significantly reduce the ionization effectiveness of the nuclear products 
are energy loss to walls and excitation without ionization. Nonetheless, for conditions typical of MHD 
concepts, namely, characteristic lengths on the order of 10-100 cm and density on the order of 10 3 — 10 
STP density, each nuclear interaction can be expected to produce hundreds to thousands of ionization 
events. As a result, a population of free electrons can form which will enhance the electrical conductivity 
through electron impact excitation and ionization. 

As with most innovative concepts, initial work can be traced back several decades. For example, the 
1 He neutron interaction was first examined during the mid- 1 960 ’s at AB Atomenergi in Sweden by Braun 
et al. 24 This experimental effort covered a temperature range of 300 to 1,600 K and density ranging from 
0.25 to 1 times standard atmospheric density with a maximum thermal neutron flux of 10 1 '/cm 2 sec. 25 
Their measurements over this limited range of conditions, however, were not promising; i.e., 0.1 S/m, and 
they concluded that the nuclear-induced conductivity enhancement was insufficient for MHD energy 
conversion applications. 
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These results discouraged further study for many years. Bitteker, however, expanding on the two- 
temperature helium plasma model of Watanabe et al., 26 has recently developed a kinetic model for Tie 
which indicates that there is a regime of thermodynamic and neutron flux conditions favorable for mean- 
ingful conductivity enhancement. 27 

Bitteker’s calculation for 7 He conductivity is shown in figure 13 as a function of relative density for 
different thermal neutron flux levels. It is clear from these plots that a maximum conductivity enhancement 
exists at some density as a function of neutron flux. At lower density, the rate of ionization is too low to 
produce sufficient ions to maintain the conductivity enhancement. At higher density, the rate of recombina- 
tion outpaces the rate of ionization while simultaneously the charge mobility decreases, resulting in lower 
conductivity. Nonetheless, these graphs show that nuclear-induced ionization processes can result in very 
significant conductivity enhancement provided a sufficient neutron flux. 



Figure 13. Computed electrical conductivity versus relative density for pure 'He. 
Solid lines correspond to a gas temperature of 1,500 K. 


Clearly, nonequilibrium conductivity enhancement requires a continuous input of ionization 
energy, in this case, neutrons. This fact would suggest that the entire generator must be located within the 
neutron field to operate. However, the rate at which the conductivity enhancement decays after the removal 
of the ionizing energy source is a function of the kinetics of the resulting plasma. 
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The ion-electron recombination process may take many routes, depending on the particular species 
involved, but usually some form of three -body recombination is dominant. In any case, these recombina- 
tion processes are a function of the electron energy. For many gases of interest, the three -body recombina- 
tion cross section decreases as the electron temperature increases above the bulk temperature, which 
decreases the rate of recombination, and hence, slows the decay of the conductivity enhancement. This 
effect becomes particularly pronounced when the electron temperature is above 0.5-1. 0 eV for a typical 
MHD flow below a few thousand degrees Kelvin. Under these conditions, the decay of the conductivity 
enhancement may be decreased by orders of magnitude. This quasi- stabilized state of ionization is 
sometimes termed “frozen inert gas plasma.” 

For conditions typical of concepts using nuclear-enhanced conductivity, the duration of useful 
conductivity enhancement may be extended to milliseconds to tenths of seconds after exiting the ionizing 
energy source without any additional energy input. This time may be further extended if auxiliary electron 
pumping systems such as radio frequency sources are employed to maintain the electron temperature 
during transit. 

These arguments combined with the recent kinetic calculations make a strong case for reinitiating 
a vigorous research program devoted to neutron-induced ionization reactions. In fact, NASA Marshall 
Space Flight Center staff are developing detailed research plans that address these critical technology 
issues on a firm experimental and theoretical basis. In the authors’ opinion, the potential payoffs for nuclear 
MHD space power in terms of design simplicity, enhanced performance, and increased reliability more 
than compensate the risks associated with a focused research effort. 
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5. MISSION CAPABILITY ASSESSMENT 


With the detailed nuclear MHD system analysis results inhand, it was of interest to examine how 
this technology could impact deep space mission capability. Because of its widespread benchmark appeal, 
the authors selected a piloted Mars mission as a baseline case, even though the technology is fundamen- 
tally suited for a much wider range of applications including interplanetary cargo transport, outer planet 
exploration, and routine commercial Earth-lunar transport. 

In this cursory study, attention was focused on a hypothetical spacecraft using a closed-loop nuclear 
MHD Brayton cycle powerplant coupled to high-power electric thrusters. In order to remain consistent 
with the previous point design analysis, baseline calculations assumed a 100-MW th , gas-cooled reactor 
driving a disk MHD generator having an enthalpy extraction ratio of 40% and an isentropic efficiency of 
70%. This implied a thermal efficiency of 40% and a thruster input power of 40 MW e . For near-term 
subsystem technology assumptions, this implied an overall propulsion system specific mass of 1.1 kg/kW e 
(see fig. 9). The delivered payload was assumed to be 100 t, unless otherwise noted. 

Using the widely recognized spacecraft trajectory optimization CHEBY-TOP code, several Mars 
trajectories were examined with respect to the 2018 Mars mission opportunities. The basis for this 
trajectory optimization algorithm is rooted in the early work of Melbourne at the NASA Jet Propulsion 
Laboratory. 28 

For all of the results reported here, the code was used to find the optimal departure date that mini- 
mized propellant usage for the outbound leg to Mars. A full round trip optimization could change the 
departure date and result in increased propellant requirements for the outbound journey, but the present 
study did not account for this effect. All trajectories were assumed to start in a 1,000-km circular Earth 
orbit (i.e., sufficiently high for safe operation of a nuclear reactor). The spacecraft was then spiraled out 
from the Earth, thrusted or coasted through heliocentric space, and spiraled into a 500-km Mars orbit. 

Figure 14 shows the initial mass in low-Earth orbit (IMLEO) required to deliver a 100-t payload to 
Mars for trip times of 120, 150, and 180 days. The power was held constant at the design point of 
40 MW e . These results show that high-power thrusters with specific impulse capabilities between 5,000 
and 8,000 sec are attractive for this mission scenario. Within this specific impulse range, it is observed that 
the payload can account for as much as half of the IMLEO. The outbound trajectories for the 180- and 
120-day transits are shown in figure 15. 
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Because the CHEBY-TOP code also had the capability of determining the optimal power for a 
given system configuration, it was of additional interest to examine those predictions for the 100-t payload. 
The predicted optimal power characteristics are shown in figure 16 along with the corresponding system 
specific mass and IMLEO. Here, the optimal power for delivering 100 1 in 120 days stalls out at 13.2 MW e , 
at a specific impulse of 2,000 sec, and climbs to 30.6 MW e , at a specific impulse of 10,000 sec. Thus, the 
40-MW e point design is substantially greater than the optimal thruster power requirement. 
Simply put, the power supply is oversized for the stated mission parameters. 




CD 


/ sp (sec) 


Figure 16. Optimal power characteristics of a nuclear MHD space 
power installation for a 100-t payload. 


The 40-MW e installation point design could be used optimally in two alternative ways: faster trip 
times and/or larger payloads. For example, if the trip time is held constant at 120 days and the payload is 
allowed to vary so that thruster power is optimized at 40 MW e , the results shown in figure 17 are obtained. 
In this case, it is possible to deliver 293 t to Mars in 120 days using a thruster having a specific impulse of 
2,000 sec; however, the IMLEO, in this case, is a monstrous 880 1. For a thruster having a specific impulse 
of 10,000 sec, the delivered payload drops to 132 t, but the IMLEO becomes much more reasonable at 
225 t. Note that any high-power thruster technology having specific impulses >6,000 sec will yield a 
payload fraction (mass of payload/IMLEO) greater than 0.56. This is a very attractive feature for 
fast 120-day transits. 
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Figure 17. Optimized payload for a 40-MW e nuclear MHD space power installation. 


While admittedly incomplete, these mission analyses indicate the enormous potential of nuclear 
MHD space power installations when coupled to high-power electric thrusters. The authors strongly 
emphasize that these mission capabilities are based on a practical systems level assessment with conserva- 
tive projections for subsystem technology development. 
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6. CONCLUSIONS 


After thoughtful examination of deep space propulsion requirements and potential technological 
solutions, the authors concluded that nuclear fission reactors coupled with MHD generators could serve as 
suitable low specific mass powerplants for high-power electric propulsion systems. These systems display 
the necessary energy density requirements while also remaining entirely within the realm of current engi- 
neering feasibility. The technological development risk for multimegawatt installations is only slightly 
greater than for turbogenerator-based powerplants, yet the potential payoff is immense — an order of 
magnitude or more decrease in system specific mass. 

Based on this line of reasoning, the authors pursued a preliminary system analysis of multimega- 
watt nuclear electric propulsion (NEP) in which a gas-cooled nuclear reactor is used to drive an MHD 
generator in a closed-loop Brayton cycle. The primary purpose of this exercise was to demonstrate that the 
necessaiy subsystem technologies either currently exist or have a basis in the existing technology base with 
clear development paths for attaining the required performance attributes. 

This was accomplished by using a thermodynamic cycle analysis in context with a thorough system 
analysis to explore the available design space. The results of this analysis indicated the appropriate operat- 
ing regime as well as the range of subsystem component performance needed to make the system viable. 

The major technological hurdles may be summarized as follows: Develop reliable, long-life, 
high-temperature fissile fuel elements; demonstrate an increase in the nonequilibrium MHD generator 
enthalpy extraction ratio from 30% to 40% and isentropic efficiency from 50% to 70%; develop efficient 
high-power (megawatt class) electric thruster technologies; and demonstrate that the entire system can 
attain a useful operating lifetime. A major advantage here is that the higher performance translates into a 
shorter lifetime requirement (megawatt class) i.e., days or months instead of years. It should be noted that 
improvements in space radiator technology would also be immensely beneficial, but adequate systems 
could be developed using the existing state of the ait. 

Based on these results, the authors conclude that MHD NEP technology offers a path for high- 
payoff, near-term realization with limited development risk. As such, the authors recommend that vigorous 
technology development steps be undertaken. 
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